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Abstract 


Most of the structural designer do not consider masonry infill walls during design process due to a 
lack of modeling guidelines in design standards and are treated as non-structural elements. In fact, the 
interaction effect between bounding frames and infill masonry is a complicated issue in nonlinearity 
of structures. The current seismic codes indirectly incorporate the nonlinear response of structure 
through linear elastic approach by considering the response reduction factor ‘R’ without comprising 
infill. In this context, this study evaluates the response reduction factor of existing engineered designed 
RC frame structures that are designed based on Indian standard codes. For this, three existing RC 
buildings were selected and performed non-linear pushover analysis. The structural response was 
examined in terms of natural period, base shear, strength, stiffness, ductility and response reduction 
factor. The results specify that the buildings with infill walls significantly influence on ‘R’ value of 
structures. Additionally, study shows that the variation of *R" value mainly depends on the percentage 
of infill inclusion. 


Keywords: Infill wall, Response reduction factor, RC building, Ductility factor, Overstrength, Non- 
linear analysis. 


Similarly, as presented in Fig. 1, the seismic 
hazard map of Nepal shows the peak horizontal 
acceleration at bed rock for Pokhara is up to 400 
gals [2]. The peak ground acceleration at rock and 
soil sites in western Nepal were found 0.29 and 
0.39g, respectively, at 10% probability of 
exceedance in 50 years [3]. 


1. Introduction: 


In the past, Nepal has experienced 7 catastrophic 
earthquakes in different epicentral locations of the 
country. More recently, the massive Gorkha 
earthquake of magnitude 7.7 (Mw) occurred on 25 
April 2015 caused around 9000 casualties, 20,000 
injuries, and more than 500,000 houses were 


destroyed (USGS 2018; Bilham 2015). In fact, 
Nepal is located in the highly seismically active 
region. The global seismic hazard map 
highlighted Nepal in a high seismic zone (zone V) 
with possible shaking of MMI IX or above with a 
10 % probability of exceedance in 50 years [1]. 


Reinforced concrete (RC) building with masonry 
infills is the most common construction practice 
all around the world. Often, engineers do not 
consider masonry infill walls in the design process 
due to a lack of modeling guidelines in design 
codes and are treated as non-structural elements. 
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In fact, the interaction effect between bounding 
frames and infill walls is a complicated issue 
when nonlinearity response occurs. It is 
appropriate if the frame and infill panels are taken 
separately by considering enough gap in-between 
[4]. Yet, this approach is rarely used in common 
practice due to difficulties in maintaining uniform 
gap between the wall and frame. In fact, masonry 
infill can drastically modify the structural 
response and make the structure more vulnerable 
if the influence of infill is not properly considered 
during design and construction [5]. However, the 
contribution of infill to the building response can 
be positive or negative, depending upon its 
configuration, position, strength and stiffness of 
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the frames and infills [6]. As indicated in Fig. 2, 
introduction of the masonry infill in RC frame 
changes the lateral load-transfer mechanisms of 
the structure from predominant frame action to the 
predominant truss-action causing reduction in 
bending moments and increase in axial forces in 
the frame member [7]. During earthquake, the 
large in-plane stiffness and strength of the infills 
wall lead to a substantial increase in the 
robustness of building which could be regarded 
globally as the advantage of these walls. On the 
other hand, the induced local phenomena such as 
short-column effect, soft-story effect, torsion, and 
out-of-plane collapse could be the disadvantage of 
these walls (see Fig. 3). 
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Figure 1: Seismic hazard map of Nepal [2]. Bedrock peak ground-horizontal acceleration 
is calculated for five hundred years to return period. Contour- interval of 100 gal 


Figure 2: Change in lateral load transfer mechanism due to masonry 
infill (a) Predominant frame action;(b) Predominant truss action [7] 


Pokhara is an important tourist city as well as the 
capital of the Gandaki province makes it under 
intense pressure of rapid urbanization. Looking at 
the urbanization scenario in terms of the building 
construction, if the constructed buildings are not 
well designed and detailed as per the design codes, 
the consequences could be more devastating and 
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vulnerable during earthquakes. In such condition, 
a detailed assessment of the existing buildings in 
Pokhara is necessary. It gives a tentative idea 
about the current construction practices and 
responses of the structures in different level of 
ground shaking. 
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Most of the current seismic design codes 
indirectly incorporate the nonlinear response of 
structure through linear elastic approach by 
considering the response reduction factor ‘R’ 
without comprising infill. This leads to structural 
response deviating radically from what is 
expected in the design because of the incorrect 
estimation of the design lateral force. The seismic 
loads imposed in the structures are normally 
higher than the loads considered during design 
process. Most of the earthquake design standards 
today include the non-linear structural response 
through response reduction factor. In this 
situation, the present study estimates the actual 
‘R? value for engineered design reinforced 
concrete buildings in Pokhara valley. The 
overstrength and ductility reduction factor of 
representative building structures are determined. 
For numerical analysis, non-linear pushover 
analysis is performed. Three representative 
engineered designed buildings with  infilled 
masonry walls are selected and analysed. Finally, 
the structural response was examined in terms of 
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Figure 3: (a) Soft story failure; (b) Short column effect due to window opening; (c) and (d) Short 
column effect due to partial height of masonry infill wall [8] 


natural period, base shear, strength, stiffness, 
ductility and response reduction factor. 


2. Methodology for Calculation of Response 
Reduction Factor: 


The response reduction factor 1s used in modern 
seismic codes to scale down the elastic response 
of a structure. This factor allows the designer to 
use linear elastic forced based design while 
accounting for nonlinear behavior and 
deformation limit. The concept of the ‘R’ factor 
is based on the fact that a well-detailed frame 
system could sustain the large inelastic 
deformation without a collapse in a severe 
earthquake and develop lateral strength in excess 
of their design strength [9]. 


The response reduction factor was first proposed 
by the ATC-19 [10]. It is described by the product 
of the three factors i.e. ductility factor, 
overstrength factor, and added viscous damping. 


R = R, X R; X Rz (1) 
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Figure 4: Base shear force vs roof lateral displacement with structure response parameter (a) 
Proposed ‘R’ by ATC-19 [10] (b) Proposed ‘R’ according to IS-1893 [12] 
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In Eq. (1) ‘Ry is the ductility factor, ‘R,’ is the 
strength factor, ‘Rọ is the any added viscous 
damping. In the absence of any supplemental 
viscous damping, the ‘Rọ is taken as unity. 
Further ATC-34 [11] accounts for the redundancy 
of the structure i.e. redundancy factor, which 1s 
mostly used in the current scenario for the 
formulation of the response reduction factor. The 
evaluation of the R- factor is based on the 
formulation proposed by the ATC-34 [11]. 


R =R, X Rs X Rp (2) 


In Eq. (2) ‘Rp is the ductility factor, ‘Rẹ is the 
strength factor, ‘R,’ is the redundancy factor. The 
assessment of this factor can be obtained from the 
pushover curve or base shear force vs roof lateral 
displacement relationship (Fig. 4) 


2.1. Ductility Factor: 


Ductility factor measures the ductility of a 
structure as a whole. It depends on the expected 
response during earthquakes which is associated 
on the design and detailing criteria. Newmark and 
Hall [9] first tried to relate the ductility factor 
with displacement ductility ‘p’ for a single degree 
of freedom system. The relation of ‘R,’ proposed 
by Newmark and hall [9] is used in the study. 


R, = 1.0 for T < 0.2 sec (3) 


R, = 42u — 1 for 0.2 < T < 0.5 sec (4) 


R, =u for T > 0.5 sec (5) 
2.2. Overstrength Factor: 


The researchers have specified that structure 
tends to have more capacity than accounted for in 
the analysis. The structure could take the forces 
considerably larger than they were designed for. 
This is due to the presence of reserve strength in 
the structure which is not accounted for in the 
design [13]. The primary sources of the 
overstrength are sequential yielding of the critical 
region, material overstrength, partial safety factor 
used for load and material, confinement of 
concrete, strain hardening, strength contribution 


of non-structural elements etc. [14]. 
Mathematically, it is expressed as. 

V 
R= Uy. (6) 


In Eq. (8) ‘Viy’ is the first yielding strength and 
‘Vg’ is the design strength. The ATC-34 [11] use 
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over strength and redundancy factors as two 
separate parameters, but some literature 
considered redundancy as a _ parameter 
contributing to over strength and assume both as a 
single unit [15]. In this study, it is assumed 
redundancy and over strength factor as a single 
unit in such condition over strength ratio 
becomes: 


V 
R= b, () 
2.3. Formulation Used in This Study: 


All the buildings were analyzed and designed 
with IS-1893 [12] which provides a real force for 
elastic structure and divides that force by 2R to 
account for the nonlinearity of the structure. To 
evaluate ductility and overstrength values from 
the  force-displacement relationship curve, 
bilinear idealization of the curve 1s needed. 


(Elastic strength demand) 


2R = SRNR, 


(Design strength) 
2.4. Modelling of Frame Element and Infill 


Panel: 
2.4.1. Beam and Column: 


The finite element program SAP 2000 [16] 
provides section designer feature to create the 
section of existing beam and column with 
reinforcement detailing but it cannot assemble all 
the sections having different reinforcement 
detailing into single unit. Since plastic hinges are 
formed at the most stressed section i.e. support 
section, we have to use support section detailing 
throughout the whole length of frame member. 


2.4.2. Infill Panel: 


To perform the safety assessment of existing RC 
infilled framed structures or to design new ones, 
the nonlinear behavior of these structural systems 
during a strong earthquake should properly be 
taken into account [17]. 


There are many different techniques proposed in 
the literature for the simulation of the infilled 
frames, which can be divided into two groups; 
namely micro-model and simplified macro 
models. The micro model considers the high level 
of discretization of the infill masonry panel, in 
which the panel is divided into numerous 
elements to take into account the local effect in 
detail. Whereas, the macro model is most 
commonly used technique due to its simplicity. In 
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this model, infill panel is modelled as a diagonal 
strut with strut width and the thickness same as 
that of the wall. Different researchers have 
proposed different formulae to calculate the strut 
width. The method proposed by the Paulay and 
Priestley [5] was used in this study where the 
strut width is taken as one fourth of the diagonal 
length of masonry wall. The diagonal strut is 
usually connected to the intersection points of the 
beam and column centerlines. This strut is 
designed in such a manner that it only carries 
compression. Transfer of bending moment from 
RC frames to masonry was prevented by 
specifying the moment release at both ends of the 
strut, the moment release is accompanied by 
using pinned joint at the ends of the strut. 


Jkl 


AN 
Aoo dg 


Figure 5: Typical infill modeling [18] 


Displacement 


Figure 6: Force-displacement relationship for the 
equivalent strut model [19] 


For performing nonlinear pushover analysis, 
user defined plastic hinges are assigned to the 
center of the diagonal strut which incorporate 
the strut behavior beyond the elastic level. For 
this, determination of the strength and the 
stiffness of the masonry beyond the elastic level 
is necessary. The constitutive relation proposed 
by Panagiotakos and Fardis [19] was used in this 
study to address the nonlinear behavior. The 
nonlinear behavior is characterized by a multi- 
linear envelope curve. 


These curves mainly have four segments. The 
first part shows the initial shear behavior of the 
un-cracked panel. The second section depends 
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on the equivalent diagonal strut formation on the 
panel after the separation of the infill from the 
bounding frame. The third part defines the infill 
wall’s softening behavior after the critical 
displacement ‘Sm’ and characterized by the 
*K3' slope. The fourth part, horizontal line, 
describes the final behavior of the infill walls. 
The final part can also be ignored by assuming a 
softening line that reaches a zero residual 
strength ‘Su’ as shown in the Fig. 6. 


3. Description of the Case Study Building: 


In this study, three existing RC-MRF buildings 
were randomly selected and designed based on 
the IS-1893 [12] considering PGA value of 0.36 
g that corresponds to zone-V of Indian seismic 
hazard map. The buildings have a maximum of 
4 stories with story height 3 m and only the 
portion above the plinth level are considered 
during the analysis. 


3.1. Material Properties and Loading: 


The staircase dead load (DL) and live load (LL) 
along with wall load (with and without opening- 
considering 30 %) is transferred to respective 
floor beam as equivalent UDL. Seismic loads 
are considered acting in the horizontal direction 
(either of the two principal directions) and not 
along the vertical direction. The floor 
diaphragms are assumed to be rigid. The unit 
weight of concrete and brick including plaster 
are assumed 25 kN/m? and 20 kN/m’, 
respectively. Poisson’s ratio for concrete and 
masonry be 0.2, dead load on floor finish = 1.2 
kN/m”, live load on floor = 2 kN/m? (25 % for 
earthquake), live load on roof = 1.5 kN/m? (nil 
for earthquake, IS 1893:2002 [12], dead load in 
staircase = 4 kN/m”, live load in staircase = 3 
kN/m’, full brick wall load = 12.36 kN/m, half 
brick wall load = 6.18 kN/m. All the dead loads 
and live loads values were taken from the IS 875 
[20, 21] respectively. The material properties of 
all three buildings are identical throughout all 
stories. Assuming the yield strength of 
reinforcing steel 500 Mpa and compressive 
strength of concrete 20 Mpa. 


3.2. Description of Structures: 


Due to the geometrical irregularity, the structure 
is analyzed with a three-dimensional (3D) 
model. The beam and column sections are taken 
as (230x430) and (300x300) mm, respectively. 
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The sectional size and reinforcement detailing of 
the plinth beam are identical for all the three 
buildings (230x300) mm and (3-166 top, 3-166 
bottom), the thickness of slab-100 mm, exterior 
and interior wall are 230 mm and 115 mm 


c 
NU LL 
Waa" 


(a) 


respectively. The plan and 3-D of case study 
building models are presented in Fig. 7. 
Similarly, Table 1 summarized beam-column 
sectional dimension and reinforcement 
detailing. 


(b) (c) 


Figure 7: Plan and 3-D model of building: (a) Type 1, (b) Type 2, (c) Type 3 


Table 1: Beam column dimension and reinforcement details 


Seis Story Member vers em Reinforcement details 
Cl 300X300 = 4-20 + 4-16 
Ist C2 300X300 8-16 
1 B 230X430 3-16(top), 2-16h+1-12(bottom) 
Cl 300X300 8-16 
Es os C2 300X300 6 -16+2-12 
B 230X430 3-16(top), 2-16h+1-12(bottom) 
All floor C 350X350 6-20+2-16 
2 (1,2,3 Bl 230X430 —- 2-1 6p + 3-1 2¢h(top), 2-16d--1-12 (bottom) 
&4) B2 230X430 3-16¢h(top),2-1 6+1-12(bottom) 
Cl 300X300 8-20 
im C2 300X300 = 4-20 + 4-16 
C3 300X300 8-16 
B 230X430 . 4-l6d(top)2-16i--1-12á(bottom) 
Ci 300X300 . 4-20 + 4-16 
2nd & C2 300X300 8-16 
3rd C3 300X300 6-16+2-12h 
B 230X430 . 4-l6d(top)2-16--1-12ó(bottom) 
Top Cl 230X430 6-16+2-12h 
(cover) B 230X430 . 4-l6d(top)2-16--1-12á(bottom) 
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4. Analysis and Interpretation of the Result: 
4.1. Natural Period: 


The comparative result of the natural period with 
and without infill consideration as shown in Fig. 
8, indicates that the inclusion of the infill 
interaction decreases the natural period of the 
structure. The decrease in period ranges from 40- 
60 % of the bare frame in both directions. The 
result is consistent with the previous study [22]. 
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Figure 8: Comparative Eigen period with and 
without infilled frame in (a) X-direction and (b) 
Y-direction 


4.2. Base Shear Capacity: 


The seismic forces are imposed higher at the 
action 


higher floor level and cantilever 
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accumulates these forces from top to bottom. The 
sum of all forces acted at the ground story are 
resisted in the opposite direction and termed as 
base shear. As indicated in the Fig. 9, the 
introduction of the infill drastically increases the 
ultimate base shear of the structure which may 
subsequently increase the design base shear due 
to decrease in the period. The increased base 
shear of the studied buildings with the inclusion 
of infilled ranged from 1.2 to 2.2 times greater 
than in the case of bare frame. The result is in the 
range to the previous study [23]. 
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Figure 9: Comparative base shear with and 


without infilled frame in (a) X-direction and (b) 
Y-direction 
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Figure 10: Type 1 capacity curve with and without considering infill in X and Y-direction 
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4.3. Response Reduction Factor Value: 


The presence of infil panels increases the 
stiffness and strength of the system. This is 
confirmed by studying the pushover curves for 
the bare frame and infill masonry structure. The 
study highlights that the average ‘R’ value of the 
study building models are in the range of (3.5-4). 
The Indian standard code [24] also mentions that 
inclusion of infill may increase or decrease the 
‘R? values which ultimately increase or decrease 
the design force. In such condition, code clearly 
specifies that the design force must be calculated 
with and without considering infill and use larger 
force out of the two conditions. Moreover, the 
other researchers also reported an increase in ‘R’ 
value due to substantial increased value of 
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overstrength [25-27]. The ‘R’ value for studied 
building models are summarized in Figs 10-12. 


4.4. Comparison 
Factor: 


of Ductility Reduction 


Due to an increase in the stiffness, natural period 
and displacement ductility of the structure 
decrease. Consequently, the ductility reduction 
factor decreases because the ductility reduction 
factor is a function of natural period and 
displacement ductility (p). The decreased 
percentage of the ductility reduction factor by the 
inclusion of the lateral stiffness of the infilled is 
about 35-45 % in the both directions. Ductility 
reduction factor of an infilled and bare frame in 
both X-direction and Y-direction in shown in Fig. 
13 (a) and (b). This value is similar to the values 
observed in previous studies [26, 27]. 
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Figure 11: Type 2 capacity curve with and without considering infill in X and Y-direction 
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Figure 12: Type 3 capacity curve with and without considering infill in X and Y-direction 
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Figure 13: Ductility reduction factor of an infilled and bare frame in (a) X-direction, (b) Y-direction 
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4.5. Comparison of Overstrength Factor: 


The previous study [7] concluded that the initially 
larger amount of lateral load attracted by the 
masonry infill but when it crosses its ultimate 
load-carrying capacity, there will be a sudden 
drop in the strength of the structure, the same 


Type of Building 


0 2 4 6 
Overstrength Factor 


pattern is seen in the Figs. 11, 12 & 13. From the 
above analysis, it is clear that unreinforced 
masonry infill frame has an average of 1.5 to 2 
time's higher strength than the bare frame. 
Overstrength factor of an infilled and bare frame 
in in both X-direction and Y-direction in shown 
in Fig. 14. 
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N 


— 


Type of Building 


0 2 4 6 8 
Overstrength Factor 


Figure 14: Overstrength factor of an infilled and bare frame in (a) X-direction (b) Y-direction 


5. Conclusions: 


The empirical expression for the determination of 
time period as mentioned in standard codes shows 
a lower value than the one calculated by the finite 
element approach. This shows that the code 
formula for estimating the time period doesn’t 
account for the irregularities present within the 
building. Hence the code formula for determining 
the time period is insufficient to include all the 
responses when subjected to the earthquake 
excitation and should be amended. The study of 
mode shape shows that for regular or nearly 
regular building, the first and second modes are 
purely translational and the third mode is purely 
rotational. While as with the increase in the 
irregularity of the building the mode shapes 
become completely of mixed type. The study 
shows that with the increase in the irregularity of 
building, the number of modes to achieve 90 % of 


modal mass participation for accurate 
determination of dynamic response increases. 
Under the unidirectional application of 


earthquake loading, the response induced in the 
perpendicular direction to earthquake excitation 
increases with increases in the plan irregularity 
due to lateral-torsional coupled behavior of the 
irregular buildin This study evaluates the response 
reduction factor of the existing RC buildings with 
masonry infills. In the analysis, the influence of 
gravity load is included but the P-A effect and 
soil-structure interaction are ignored during 
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numerical analysis. The results clearly indicate a 
significant variation of fundamental time period, 
base shear, ductility, overstrength, and response 
reduction factor in the presence of infill masonry 
wall. Thus, considering the infill panel in the 
analysis would influence the seismic behavior of 
frame structure to a greater extent. Consequently, 
‘R’ factor either decreases or increases depending 
upon the relative increase or decreased values of 
the overstrength and ductility factor. The main 
conclusions are summarized as: 


= the presence of infills decreased the time 
period of structure by 40-60 % compared 
to the bare frame. 


= (the inclusion of infills masonry increases 
the base shear from 1.2 to 2.2 times than 
bare frame. 


= introduction of infilled decreases the 
displacement ductility factor by 35-45 96. 


= for infill opening greater than 30 96, 
decreased percentage of overstrength was 
relatively higher than the increased 
percentage of displacement ductility, 
resulting into a decreased value of the 
response reduction factor. 


= generally, the infill opening is considered 
in the range of the 30 % and above for 
residential building. 
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for the design of residential building, the 
study recommend the ‘R’ value in the 
range of 3.5-4. 
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